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We study the critical behavior of lattice Quantum Chromodynamics (QCD) in the strong 
coupling approximation with Kogut-Susskind and Wilson fcrmions at finite temperature 
(T) and zero chemical potential. Using the Hamiltonian formulation we construct a mean 
field solution to the equation of motion at finite T and use it to study the elementary 
thermal excitations and to extract some critical exponents characterizing the observed 
second order phase transition. We find similar critical behaviors for Kogut-Susskind and 
Wilson fermions at finite T. 
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In this letter we study the critical behavior of lattice Quantum Chromodynam- 
ics (QCD) in the strong coupling approximation"'^ using Kogut-Susskind and Wilson 
fermions at finite temperature (T) and zero chemical potential. This work is a nat- 
ural continuation of Ref. El where the strongly coupled lattice QCD was examined 
at zero T but at finite chemical potential. We shall follow the same method as in 
Ref. 121 and use the Hamiltonian formulatiorP to construct a mean field solution to 
the equation of motion at finite T . This solution is then used to study the elemen- 
tary thermal excitations and to evaluate some critical exponents characterizing the 
observed second order phase transition. Comparisons of critical exponents as well 
as of thermodynamic quantities indicate that Kogut-Susskind and Wilson fermions 
behave similarly in the neighborhood of their critical points at finite T . 

Analytical studies of strongly coupled lattice QCD at non zero temperature (and 
chemical potential) have been presented previously by several authors both in La- 
grangian and Hamiltonian formulations using mostly Kogut-Susskind fermions.El 
Where the behavior of the theory has been investigated for finite T and zero chem- 
ical potential, a second order phase transition has been observed with the chiral 
condensate playing the role of the order parameter. The order parameter is fi- 
nite for temperatures lower than some critical temperature {Tc) and vanishes for 
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T > Tc- This is because the Lagrangians and Hamihonians that were studied have 
the structure of an anti-ferromagnet with a nearest neighbor interaction, and for 
zero current quark mass this system is known to exhibit a second order phase tran- 
sition at finite In this letter we examine more closely the critical behavior of 
the observed second order phase transition using both Kogut-Susskind and Wilson 
fermions. 

As in Ref.lSlour starting point is Smit's effective lattice Hamiltonian for strongly 
coupled QCD consisting of only quark fields.^ This Hamiltonian has been derived 
using second order strong coupling perturbation theory and only involves terms 
to order 0{l/g'^) and 0(1/-/Vc), where g and Nc are bare quark-gluon coupling 
constant and the number of colors, respectively. These terms are bilinear in quark 
fields. Terms which are trilinear in quark fields appear at order 0{l/g'^) for TVc = 3 
although they are not included in the Hamiltonian of Ref. El We stress that our 
method presented here is independent of the structure of the Hamiltonian, and 
is therefore applicable even in the presence of these 0{1/ g'^) terms. The effective 
Hamiltonian is expressed in the temporal gauge supplemented by the Gauss's law 
constraint. This constraint restricts us to work in the subspace of the lattice (single 
site) Hilbert space where all physical states in the space are SU(3) color singlets 
and therefore locally gauge invariant. 

Adopting the notation of Smit, El setting the lattice spacing to unity and working 
in momentum space, the charge conjugation symmetric form of Smit's effective 
strongly coupled ^ 1) QCD Hamiltonian is 



K 
"SiVc 



«((P2+P3)-"l) _|_ g«((Pl+P4)-n!)l 



pi,--,P4 ; 
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(1) 



where (S;) = — i (707; — *f7o) with the Wilson parameter r taking on values between 
and 1. The quark field is a Kogut-Susskind fermion for r = and a Wilson 
fermion for < r < 1. Using the standard summation convention color, flavor 
and Dirac indices for the fermion field are denoted by (a, 6), {a, (3) and (p, 1^,7, (5), 
respectively. Mq is the current quark mass and the coupling constant K is related 
io g as, K = 2Nc/ (N^ — 1) l/g^. Therefore model parameters of i?off are the Wilson 
parameter r, the current quark mass Mq and the effective coupling constant K. 

For Kogut-Susskind fermions and for vanishing current quark mass (r — Mq = 
0) the effective Hamiltonian possesses a U{4Nf) (global) symmetry where Nf is the 
number of flavors. Smit has found that this symmetry is spontaneously broken down 
to U{2Nf) (E) U{2Nf) accompanied by the appearance of 8Nj Goldstone bosonsPl 
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A finite current quark mass also breaks the original U{ANf) symmetry down to 
U{2Nf) (g) U{2Nf) albeit explicitly. For Wilson fermions the latter symmetry is 
further explicitly broken to U{Nf) thereby solving the fermion doubling problem. 

The efi^ective Hamiltonian provides a description of only the ground state of 
strongly coupled lattice QCD. In this state there are no excited color electric flux 
links. Therefore the external parameter T must be small enough to not to thermally 
excite any color electric fluxes. This condition is satisfied for T < 1/kW 

In Ref.0we presented an ansatz satisfying the equation of motion corresponding 
to Hcff for all values of T and fi and explored the consequences in the T — s- limit. 
In this work we shall take the other limit of /i ^ 0. In this limit the ansatz takes 
on a simple form given by 



(2) 



In Eq. Q the thermal field operators B and annihilates a particle and creates 
a hole, respectively, while D and are the annihilation and creation operators for 
an anti-particle and an anti-hole, respectively. These operators satisfy the free field 
fermion anti-commutation relations and each annihilation operator annihilates the 
interacting thermal vacuum state \Q{T)). The appearance of holes and anti-holes 
reflects the doubling of the Hilbert space inherent in the formalism of thermo field 
dynamics which we have used in Ref . 

The coefficients Up and /?p are = (1 — UpY^'^ and f3p — (rip)^/^ where Up is 
the Fermi distribution function Up = [c'^p/C'bT) _|_ ^j-i -v^jth ujp — sin^(p • hi) + 
M{p )^Y/^ . As will be explained shortly M(p) is to be interpreted as a dynamical 
mass. It is in general both momentum and temperature dependent and is the only 
unknown quantity, and not a parameter, in our ansatz of Eq. Finally, spinors 
^ and rj satisfy the equation of motion of a free massive lattice Dirac Hamiltonian 
i^o given by 



^lit,p),^,it,p) 



(3) 



The main idea behind our approach is to use the equation of motion for our 
ansatz which may be expressed as 



{^ao.)^{t,p),H' 



(*aa)^(t,p),i?cff 



(4) 



to derive a self-consistency equation satisfied by the unknown quantity M{p ). Here 
the symbol : : denotes normal ordering with respect to the thermal vacuum \Q{T)). 
Our solution for the strong coupling theory defined by Smit's Hamiltonian is ob- 
tained by solving the self-consistentcy equations for M(p) for each value of T. 
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With our ansatz of Eq. |5J obeying the free massive lattice Dirac equation, Eq. @ 
describes the quark field in as a traditional quasi-particle with an effective, 
or dynamical, mass M(p). Hence the operators B\ B\ and D'^ are to be in- 
terpreted as creation operators for quasi-particles, quasi-holes, quasi-anti-particles 
and quasi-anti-holes, respectively. 

After normal ordering with the aid of computer algebra, we obtain the following 
equation of motion from Eq. Q 

fVsin(p-nO(7o7/)+^^(p)(7o)l {■^ao)sii.v) = 

I 

[^A,(p)(7o7z) + S(p)(7o)l (vi/„o),(i,p). (5) 

Thus the dynamical mass M(p ) is given by the coefficient B{p ) which is found to 
be 

B{p) - Af (p) 



I K ^ , M(q) 



2N, 
q 



[l — 15r^) cos(q • hi) cos(p • hi) 



(l + r^) sin(q • hi) sin(p • hi) 



(6) 



where the chiral condensate (^'5} is given by 

imk " 5F^<«ml[(*")-(*-.)"l-l»m> 

p '^P 

We find that the ansatz of Eq. Q and the equation of motion allow us to express 
the temperature dependent part of the vacuum energy density in the following 
deceptively simple form 



1 



2VNfNc 



i?eff)T>0-(i/eff)T=oJ = Y^^^^Mo + iu{p) ^ (3^ (8) 

According to this definition the energy density is positive and does not depend 
explicitly on the strength of the interaction K. However we should recall that the 
dynamical mass that enters in Eq. (jHl not only depends on K but also on r and 
T. The expression Eq. JH)) forms the basis of the thermodynamics of the strongly 
coupled theory under investigation. 

We also find the off-diagonal effective Hamiltonian to second order in thermal 
field operators to be 

(i?e4ff|OT) = 
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J2 [c^pPp ■Hi (p )bL.(-p )sL(p ) + c^pPp n2 (P (p )dU-p ) 
p 

+ al Hsip )^L(-P )^L(P ) + «p/3p W4(P )dUp )bU-P )] I GiT)) (9) 

where 7ii (p ) are complicated functions to be presented elsewhere. In Eq. the 
elementary thermal excitations consist of particle-hole (B^B^), anti-particle-anti- 
hole particle-anti-particle (D^i^t) a^d hole-anti-hole {D'^B'') excitations. 
Note that these excitations, which carry zero total three-momenta, are color singlets 
and therefore satisfy the Gauss's law constraint imposed on iJcff- 

In the strong coupling limit the ground state of the QCD lattice Hamiltonian is 
any state with zero gauge flux implying that it is infinitely degenerate. Smit's ef- 
fective Hamiltonian is obtained by applying perturbation theory about this ground 
state. The lowest nontrivial contribution comes from second order perturbation 
theory and correspond to propagations of mesonic objects on the lattice. In our for- 
mulation these propagations are represented by the elementary thermal excitations 
in Eq. Q and are responsible for lifting the infinite degeneracy of i?eff ■ 

Thus far the results obtained using our ansatz Eq. (0) have been exact. We shall 
now invoke the mean field approximation where the equation for Af(p), Eq. ©, is 
solved only to 0{N^) resulting in a self-consistency equation with a structure simi- 
lar to that of the Nambu-Jona-Lasinio model in the Hartree-Fock approximation.^! 
As a consequence, the dynamical quark mass, and therefore the chiral condensate 
through Eq. j?)), become momentum independent. We note that in all of the pre- 
vious analytical studies of strongly coupled lattice QCD at finite equivalent 
approximations have been invoked resulting in momentum independent chiral con- 
densates. 

We shall also work in the limit of vanishing current quark mass Mq = 0. In Ref.El 
Smit bosonized the effective Hamiltonian, identified the quark bilinears as funda- 
mental meson fields and showed that for Wilson fermions the pion mass vanishes 
when Mq — 6Kr^ . This " critical" value of the current quark mass was then defined 
to be the chiral limit of the theory with Wilson fermions. Consequently, however, 
Le Yaouanc et.al. using the same Hamiltonian as Smit claimed that above critical 
value of Mq results in a vacuum which is not chirally degenerate, leading to the 
conclusion that the massless "pions" obtained by Smit are not Goldstone bosonsP 

The method presented here avoids the bosonization method and explicitly con- 
structs a solution to the equation of motion using only fermionic degrees of freedom. 
Furthermore, our ansatz is non-confining since it is a solution of a free massive lat- 
tice Dirac equation. Thus it does not make sense to determine the pion mass, nor 
masses of any other bound states, using this ansatz. However the ansatz is not un- 
reasonable since it exactly diagonaliz es th e effective Hamiltonian to second order in 
field operators when T = and > O.EEIln addition, our formalism can accommo- 
date a description of bound states using an improved ansatz as discussed in Ref. |2 
We shall not dwell on Smit's critical value of the current quark mass for Wilson 
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Fig. 1. -(>I"I')/(2ArfiVc) vs. T for Kogut-Susskind (r = 0) and Wilson (r = 0.25) fermions with 
vanishing current quark mass Mq = and for various values of effective coupling constant K. 
Critical temperatures are shown in Table [Tl 



Table 1. Critical temperature Tc and coefficients a, a' , /B and rj 
for various values of input parameters. Values of g corresponding 
to each K are also shown for Nc = 3. 



r 


K 


9 


Tc 


a 


a' 


P 


V 


0.00 


0.80 


0.97 


0.295 


0.023 


0.0003 


0.515 


f5i 0.000 


0.00 


0.90 


0.91 


0.438 


0.363 


0.011 


0.517 


0.000 


0.25 


0.90 


0.91 


0.364 


0.143 


0.005 


0.508 


Si 0.000 


0.25 


1.00 


0.87 


0.483 


0.596 


0.026 


0.541 


Si 0.000 



fermions until this improvement has been implemented. 

We make the standard choice of letting the chiral condensate Eq. 10 be the order 
parameter of the theory. In Fig. ^we show the behavior of ( ) as a function of 
temperature for various sets of parameters. We observe a similar quahtative behavior 
between Kogut-Susskind and Wilson fermions where the chiral condensate, being 
finite at T = 0, vanishes smoothly at some critical temperature Tc characteristic 
of a second order phase transition. For A/q = 0, the effective Hamiltonian with 
Wilson fermions possesses a U{Nf) symmetry for T < Tc and U{4:Nf ) symmetry 
for T > Tc- Table ^ shows that for each type of fermion, Tc increases with the 
coupling constant K confirming the observations made in Ref. I* In addition the 
critical temperatures shown in the table satisfy the condition Tc < l/K so that no 
electric fluxes are thermally excited. 

The behavior of the chiral condensate near Tc is best characterized by the 
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critical coefficient f3 which in this work is defined as 

y-l=Ap{Tc-rf, Tc>T (10) 

where is a constant. As shown in Tabled the critical coefficients /? are approx- 
imately equal with values close to the mean field result of /3 = 0.5. We may also 
study the response of the chiral condensate to the changes in the dynamical mass, 
and define a critical exponent t] describing this response near T k, Tq as 

S^^/SM^A/Tc-T)-" Tc>T (11) 

with being a constant. For all cases we find a linear relationship between ( ^"I' ) 
and M leading to a critical coefficient of ry w 0. 

As mentioned above, we may obtain all the thermodynamical information from 
the vacuum energy density shown in Eq. ((SJ). The qualitative behavior of the vacuum 
energy density for all of our model parameters are found to be the same. As an 
example typical behaviors for the specific heat at constant pressure Cp for both 
types of fermions are shown in Fig. 12 The results were obtained as follows. First 
the dynamical mass is determined for discrete values of T by solving Eq. © to 
0{N^). In solving Eq. ^ the mode sum is replaced by a contimmm integral as 
has been done in Ref. |3 The resulting dynamical mass is then used to evaluate the 
vacuum energy density using Eq. ^ followed by a polynomial fit to approximate 
the energy density. Finally, the interpolated energy density is used to determine the 
specific heat by numerically evaluating its second derivative. 

We find the specific heats to be discontinuous at Tc- Their critical behaviors on 
both sides of the critical point are described by critical exponents a and a' defined 
as 

_jA^{Tc-Tr- Tc>T 

Ua'(r-rc)-"'T>rc ^ ^ 

with constants Aa and Aa' . The results for a and a' presented in Table Q vary 
widely, but at least they are all positive and increase with K for each type of 
fermion. 

To summarize, we began with a well defined Hamiltonian describing the ground 
state of strongly coupled QCD and introduced an ansatz from which we constructed 
a solution satisfying the equation of motion for finite T. Using the chiral condensate 
as the order parameter we find a second order phase transition for both Kogut- 
Susskind and Wilson fermions. In all cases we find that the critical temperature Tc 
obeys the condition Tc < l/K so that no electric fluxes are excited near the critical 
point. Our ansatz allowed us to study the second order off-diagonal Hamiltonian 
which clearly showed the propagation of thermally excited color singlet mesonic 
objects responsible for lifting the infinite degeneracy of the Hamiltonian. In addition, 
we have extracted some critical exponents to study the critical behavior near Tc- 
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Fig. 2. Specific heat at constant pressure Cp vs temperature T for Kogut-Susskind {r=0.0, 
K=0.80) and Wilson (r=0.25, K=0.9) fermions. 

Our main conclusion is that both Kogut-Susskind and Wilson fermions have similar 
critical behaviors at finite T. 

Finally, as stated previously we are fully aware that our mean field ansatz Eq. (0) 
is nonconfining and needs to be improved. The need for an improved ansatz is sup- 
ported by a recent numerical work by Chandrasekharan and StrouthoJ^ of strongly 
coupled lattice QCD at finite T in which mean field theory predictions were found 
to fail in the immediate vicinity of the critical temperature. How to construct this 
improved ansatz, which goes beyond the mean field approximation and allows for 
the introduction of bound states, has been described in Ref. 2 and shall not be 
repeated here. However, the method presented in this letter and in Ref. 2 to intro- 
duce temperature and chemical potential is applicable regardless of the choice of 
the Hamiltonian and the ansatz. In particular it offers an economical way to extend 
the current study to the entire T — ^ plane where a tricritical point is expected to 
be encountered.!^ 
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